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Abstract 

This article discusses entanglement between two subsystems, one with 
discrete degrees of freedom and the other with continuous degrees of free- 
dom. The overlap integral between continuous variable wave functions 
emerges as an important parameter to characterize this kind entangle- 
ment. "Beam- like" entanglement and "shape-like" entanglement are con- 
trasted. One example of this kind of entanglement is between between 
the spin degrees of freedom and the momentum degrees of freedom for a 
non-relativistic particle. This intraparticle entanglement is Galilean in- 
variant. 
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1 Introduction 

The article considers the properties of entanglement between two subsystems, 
one of which has discrete degrees of freedom (DOF) and the other has continuous 
DOF. This topic therefore lies somewhere between the more well-established 
field of entanglement in discrete systems [1] and the growing field of continuous 
variable (CV) entanglement [2J. The subsystems under consideration need not 
be "systems" in the sense of two different quantum objects, like electrons or 
atoms or fields. Following recent work [21 [H [5], by the term "subsystem" all 
that is implied is that there is an algebra of observables that can be divided into 
two commuting subalgebras and this can be used to specify a tensor product 
structure. 

Two recent applications of continuous-discrete (CD) entanglement are in 
analyses of the Stern-Gerlach experiment [B] and nuclear decays [7J. In Rj], 
the entanglement discussed is between the spin DOF and momentum DOF 
of a single particle. In [7J, the entanglement is between the total momentum 
DOF and the total spin DOF of two particles. One goal of this work is to 



put these analyses into the broader context of CD entanglement and identify 
commonalities. 

First, we will consider this problem mathematically without reference to 
any particular physical system. In this context, the importance of the overlap 
integral between the different CV wave functions will become apparent. Gen- 
erally the smaller the overlap, the greater the entanglement. There are two 
fundamentally different ways that small overlap can occur: (1) The wave func- 
tions could have support in different domains of the continuous variables. We 
use the description "beam-like" for this kind of CD entanglement. (2) The wave 
functions could have support in similar domains of the continuous variables, 
but the details of their structure lead to cancellations. We call this kind of CD 
entanglement "shape-like" . In the case of zero overlap, the different CV wave 
functions act just like discrete DOF. 

As an example, we will consider a form of CD entanglement that can arise 
in particles: spin-momentum entanglement of a non-relativistic electron. This 
kind of intraparticle entanglement is invariant under Galilean transformations 
between inertial reference frames. Therefore, it has an intrinsic physical mean- 
ing. In the concluding example, we will exemplify the two different extremes 
of CD entanglement considering the properties of Gaussian momentum wave 
functions. 

2 Continuous-Discrete Entanglement 

We will consider a system that has a complete set of commuting observables 
(CSCO) that partitions into three sets. One set will be the observables 

{P 1 ,P 2 ,...,P d }. 

These observables will be represented as operators with a continuous spectrum 
in the system's Hilbert space representation. Examples are the momentum 
generators of the Heisenberg group in rf-dimensions or the generators of trans- 
lations in R d . Another "set" will be the single operator Si, which will have a 
discrete spectrum in the Hilbert space representation. A simple example of Si is 
a component of the angular momentum, but this notion could be generalized to 
other compact Lie algebras or to a set of commuting observables with discrete 
spectra. Finally, there are the generalized Casimir observables. These will be 
represented as multiples of the unit operator in the system's Hilbert space and 
although they typically are important for characterizing the physical system, 
they will play no role in this analysis. 

From the existence of this kind of CSCO, it is implied that the system's 
Hilbert space can be factored as H — H s ®H p . To the discrete DOFs corresponds 
the n-dimensional Hilbert space H s — C". The Hilbert space associated to the 
d-dimcnsional continuous DOFs is the space of Lebesgue-integrablc functions 
Hp = L 2 (R d ), although in practice, we will limit ourselves to a dense "well- 
behaved" subspace S p c L 2 (R d ) (see below). 

We can define a basis |p,x) — |x) ® |p) f° r W in terms of a discrete label 
X € {0, 1, n — 1} and continuous label p £ R d . These are the eigenvectors of 
the operators {Pi, ...,P<j} = P and Sf. 

P|p,x) = p|p,X>, &|p,x) = ((n-l)/2- X )|p,X> (1) 

The eigenvalues of the Si operator are written as ((n — l)/2 — x) for later 
application when it will be interpreted as spin component along some axis. 

Because P is an unbounded operator, the basis eigenkets |p) are not elements 
of Hp and can only defined as anti-linear functionals on nuclear subspaces of 
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H p [S]. One suitable choice for this subspace is the Schwartz space S p of well- 
behaved vectors, i.e. vectors corresponding to smooth, rapidly-decreasing wave 
functions. For vectors in a nuclear subspace \<j>) G H. s ® S p , the basis eigenket 
expansion holds: 

n . 

10) =5Z d d P<P x (p)\P,X), (2) 
where the eigenket normalization is 

(P,X|P',X')=^X'^(P-P')- (3) 

Because of the restriction to the Schwartz space and the normalization, the 
integration measure d d p can be chosen as the standard Riemann integration 
measure in d-dimensions [5], e.g. d d p = Uf =1 dpi. 

As a consequence of the product structure of the Hilbert space, one can talk 
about internal entanglement between the discrete and continuous DOFs. For a 
pure state |</>), one can define the reduced discrete density matrix p s by 

p s = Tr p (|0)<^|) 

= E(/ rff V x (p)^(p))lx)(x'l 

x,x' 

= E^x'lxXx'l- ( 4 ) 
x,x' 

where h XfX t are the integrals of overlap of CV wave functions and the matrix 
elements of the reduced discrete density matrix. Also, one can define the reduced 
continuous density matrix (operator) p p by 

p p = Tr s (j^><0|) 

= Jd d pd d p'^4, x (p)4> x (p')^\p)(p'\ 

= |dW/(p,p')|P>(p'l- (5) 
Assuming the initial pure state is normalized, 

E / rf d Pl0x(p)| 2 = l, (6) 
x J 

then both reduced density matrices are normalized 

Tr p p s = Tt sP p = 1. (7) 

Additionally, for pure states \<fi)((j)\, but not generally for mixed states, we have 

T rp g(p s ) = Tr s .g(pP), (8) 

where g is a continuous, analytic function on the domain of Hcrmitian operators 
on both Tt s and S p . This is a consequence of the fact that the matrices p s and 
p p have the same spectrum of non-zero eigenvalues [9]. 

Because finite matrices are easier to work with, it will be easier to consider 
the reduced discrete density matrix p s for calculating entanglement properties 
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of pure states. From normalization, the sum of the diagonal elements of the 
reduced discrete density matrix is unity: 



E /l x,x = l' (9) 



The off-diagonal elements h x y = h*, x are the integral of overlap between 
different wave function components </> x and <f> x > 

= f dVx(p)<&(p)- (10) 



n x,x' 

The magnitude of these off-diagonal elements are bounded from above by l/\/2 
(for any n) by combination the Cauchy-Schwartz inequality 

0<K, x '| 2 <IVxllV,x'l (11) 
and ©, which implies that 

|/*x,xllV,x'l<V2- (12) 

The amount overlap between the wave function components determines how 
entangled the state is. To see how this works, consider the extreme cases of 
maximum overlap and no overlap. 

• Maximum overlap: The first limiting case is when all components of the 
wave function have the same momentum dependence. In this case the 
wave function has the form 

<j> x {p) = c x d>{p) (13) 

where 

\c x \ 2 = 1 (W) 



x 



and 

' d d p|<Kp)| 2 = 1 (15) 



The overlap will be truly maximal in the case where c x = w x /\/n, where 
\uj x \ — 1. In any case where (fl"3|) is valid, the discrete and continuous DOF 
can be separated and the discrete reduced density matrix elements are 

h x,x'= c x c * x '- ( 16 ) 

This is the density matrix of a pure state and the entanglement as mea- 
sured by any entanglement monotone should be zero. 

Zero overlap: The other extreme is when h XtX ' = unless x — x' ■ This in 
turn breaks into two cases: (1) All <^ x (p) are identically zero except for 
one. This state is clearly separable and has no CD entanglement. (2) The 
more interesting case is when <fr x (p) and <fi x (p) are orthogonal for all XiX' ■ 
Then the continuous DOFs can effectively be treated as discrete DOFs; 
and x then plays a dual role labeling the orthogonal CV mode and the 
discrete index. 

When h xx > — 0, the discrete density matrix p s is already in diagonalized 
form, and so the entanglement depends only on the relative magnitudes 
of the diagonal elements h x . x . Maximal entanglement occurs when when 
p s is maximally mixed and therefore I^Xixl = V n - 
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In the second case of orthogonal wave functions, it may be important to 
contrast two different kinds of orthogonality. One kind, which we metaphorically 
call beam-like, is when the wave functions are orthogonal because they have 
support on disjoint regions of M. d . For example, imagine d-dimensional Gaussian 
wave functions whose central maxima are distant compared to their widths, like 
focused beams going in different directions. Hyperentangled states [10] and 
the interferometry states used in [11] are examples of these. Alternatively, one 
could imagine that the wave functions have the similar domains but are close 
to orthogonality because of cancellations in the integration, like the various 
families of orthogonal polynomials. We call this shape- like CD entanglement. 
For example, the <fi x (p) wave functions could be energy eigenstates of some 
potential well expressed in the momentum basis. 



3 Intraparticle Spin-Momentum Entanglement 

We now consider the non-relativistic particle as an example of CD entanglement. 
For non-relativistic particles, the state space of a particle can be factored into 
the direct product of two Hilbert spaces Ti e = Tt s <£> T~C P , where Tt s is associated 
to the spin DOF and H p to the momentum DOF. This separation is allowed by 

the CSCO {P, Si] M, W, S 2 }. The eigenvalues of the invariant operators for a 

single free particle {M, W, S 2 } are the mass m, internal energy W = E—p 2 /2m, 
and intrinsic spin extracted from s(s + 1). These invariants characterize the 
representation space Ti{m, W, s), and since we will consider the electron, s = 
1/2 and m = m e is a known constants, and without loss of generality we can 
consider the case W — 0. All this can be proved by explicitly constructing the 
unitary (projective) representation spaces of the (centrally-extended, universal 
covering of the) Galilean group Q of symmetries |12j . This representation space 
corresponds precisely to the solution of the two component, non-relativistic 
Schrodingcr equation for a free electron. 

As described above, we will consider the dense subspace Tti/ 2 ® S p of the 
electron Hilbert space TC e — H(m e , 0, 1/2). Then the momentum eigenkets will 
be well-defined and we will have the usual expansion theorems: every "well- 
behaved" pure state vector in |0) 6 Ti s <g) S and can be expanded as 

\<P) = E / d3 p<Mp)lp,x> 

x J 

= />p]>> x (p)|p,x>. (17) 
J x 

We now interpret P as a momentum operator and Si as the spin component 
in some direction. By choosing the Schwartz space S p for the dense subspace, 
we guarantee that the expectation values of all positive powers of momentum 
are bounded, as well as other properties assuring that the momentum wave 
functions are "well-behaved" [8j. This restriction makes physical sense and the 
importance of energy boundedness for entanglement measures of CV systems 
has been previously emphasized |13j . 

Note that the order of the integral and sum is not relevant in (fT7]) because 
the Hilbert space of the particle factors into independent DOFs. This is not the 
case for relativistic par ticle s because for them the spin component definition is 
momentum dependent [14] . 

Because of the product structure of the nonrelativistic, single particle Hilbert 
space, one can talk about the spin-momentum entanglement of a particle in a 
frame-independent way. All shifts of reference frame are Galilean transforma- 
tions and these are represented as a tensor product of unitary operators acting 
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locally on each factor Hilbert space. Here we show the explicit construction of 
the unitary irreducible representations of the Galilean group in the representa- 
tion that is implied by our choice of CSCO: 

A general element g £ Q is denoted by ten parameters: g — (6, a, v, R), where 
b is the amount of time translation, a is the three- vector of space translation, v 
is the three- vector of velocity boost, and R € SO (3) is a rotation. The projective 
representation acting on the basis for the physical subspace of TL e can be chosen 
as 

U(g)\ P ,x) = e^-ea.v-a. v+™)J2 Dl/2 (n)xx'\p',x'), (18) 

x' 

where p' = Rp + mv and D 1 ^ 2 (R) is a 2 x 2 unitary matrix representing R in 
H-i/2 (some details of the SU(2)-SO(3) little group homomorphism have been 
swept under the table here). Notice, as promised, this representation factor- 
izes into unitary operators on the discrete Hilbert space, i.e. the unitary matrix 

D 1 ^ 2 (R), and the continuous Hilbert space, i.e. the phase factor. Since Galilean 
transformations are tensor products of local unitaries with respect to the spin- 
momentum tensor product structure, any reasonable measure of entanglement 
will be invariant across different inertial reference frames. Note, however, that 
in the relativistic case, the spin-momentum tensor product structure is not in- 
variant under Poincare sym metry, as a result, "mixing" between different kinds 
of entanglement can occur |14[ I15j (although sometimes authors have not rec- 
ognized the true origin of this phenomenon). 

The reduced spin density matrix p s of the pure electron state \<p) will have 
the following form 

s ( ^0,0 ^0,1 \ / in \ 

The eigenvalues X± of this matrix characterize the entanglement properties: 



X± = -±d--(h , Q h 1 , 1 -\h , 1 \ 2 ). (20) 

For maximum entanglement, X± — 1/2. From Ijlip . the quantity hn n/ti i — |/tn i 1 2 
is always positive and less than 1/4. The upper bound corresponds to maximum 
entanglement and the lower bound to separability. 

Separability will occur in two cases: (1) Either (po(p) or (pi(p) is identically 
zero. The state is obviously separable then because the electron is polarized in 
the direction of one spin component. (2) In the other case, (j>o{p) — w^i(p) 
where \to\ = 1. Then /io,o^i,i = |^o.i| 2 and the state can be put in the form 
(fl~3"|) where w is absorbed in cq. All other \<p) will be entangled, and this entan- 
glement increases as </>o(p) and </>i(p) come closer to equal norm of 1/V2 and 
to orthogonality. 



4 Concluding Example 

Let us now consider a specific but flexible model for considering beam-like entan- 
glement. Each wave function component 4> x {p) wm t> e taken to be a Gaussian 
in 3D momentum space with central momentum k x and width a x : 

«P) = ^374 -p(-^ P - k ^)- (21) 



G 



The constant c x gives the normalization of each Gaussian wave function, i.e. 

^p|0 x (p)| 2 = |c x | 2 (22) 



and total normalization implies |co| 2 + \c\ | 2 = 1. These kinds of states could be 
useful approximations for wave packets prepared under some dynamics that have 
an explicit asymmetry with respect to a preferred projection of spin component, 
the classic example being the Stern-Gerlach device [?]• 

The reduced spin density matrix for the non-relativistic electron with wave 
functions of the form (|2~Tj) is 

P s =( l C ° l \ (23) 
1 c* c lX * |ci|- I 



where x is the Gaussian overlap function 

x = J d 3 P o (pM(p) 



2^1 ^ 2 / -2g 2 \ 

exp l^T^J (24) 

and q = |ko — ki| is the magnitude of the difference of the central momenta. 
The eigenvalues of this matrix are 



A± = i±^-M 2 | C i| 2 (l-M 2 ). (25) 



The first observation from (|25|) is that maximum intraparticle entanglement 
can only occur if the state has a spin component expectation value (Si) = 0. 
For such states, the quantity |co| 2 |ci| 2 takes its maximal value 1/4. This feature 
is a consequence of the intrinsic preferred direction in the choice of states ([2~Tj) . 
and is not a general feature of intraparticle entanglement. 

The second observation is that intraparticle entanglement increases when 
|x| — > 0, i.e. when the Gaussian wave packets have minimal overlap. This can 
be seen most easily in the case where |co| 2 |ci| 2 = 1/4; then (f2"5)l simplifies to 

\± = \±\\x\. (26) 

The quantity |x| will go to zero if either or both 2q 2 ^> a 2 + <j\ and <jq/<ji ^ 1 
(or equivalently oq/<j\ <C 1). Physically, these limits correspond to narrow wave 
packets with different central momenta or to wave packets with same central 
momentum and very different widths, respectively. Both cases can best be 
thought of as examples of beam- like entanglement. 

Beam-like entanglement is manifest in particle experiments with beam split- 
ters, e.g. the Stern-Gerlach experiment, and in quantum optics with photon. 
It is unclear whether it is experimentally feasible to produce shape-like entan- 
glement in a free, beam-like system, but perhaps the control of wave functions 
possible within quantum optics could allow such a possibility. In bound elec- 
tron systems, shape-like entanglement between vibrat iona l energy levels and 
spin states has already be achieved with trapped ions 16 1. The orthogonality 
of the bound states allows the CV DOF to be treated like discrete DOF. 
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